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Abstract 

Using the formalism developed in a previous paper we analyze the cosmo- 
logical implications of our conclusions concerning the scaling behaviour of a 
network of cosmic strings, in particular the idea that once gravitational back- 
reaction becomes important the transient scaling regime so far explored by 
numerical simulations will be replaced by a new 'full scaling' regime, with 
slightly different parameters. This has consequences for the normalization 
of the string tension on all scales. We show how the new scaling solutions 
affect the gravitational wave bound from nucleosynthesis and the milli-second 
pulsar, and also describe the consequences from large scale structure up to 
COBE scales. 
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Cosmic strings formed at an early phase transition may have interesting cosmological effects. 
To test this idea, we need to know how a network of strings evolves. In a recent publication 

we developed a detailed analytic formalism aimed at improving our understanding of 
the evolution of such a network in a /c = Friedmann-Roberston- Walker (FRW) universe. 
The approach is based on calculating the evolution of the probability, p[r(Z)], of a random 
segment of left-moving string of length / having an extension r at time t. One of our 
most important conclusions was that the regime so far studied in numerical simulations is a 
'transient scaling' regime, which will be succeeded once gravitational back-reaction becomes 
important by 'full scaling', characterized by rather different parameters. In this letter we 
address some of the cosmological consequences arising from the new results. 

In ^ we introduced three length scales ^(t),^(t) and ((t) related to the long string 
density, the persistence length along the long strings and the small-scale structure on the 
string network respectively. In terms of the mean square extension, K{1) = r^, and ( are 
defined by 

K{i,t) ^2^{t)i, i:>t, (1) 

whereas ^(t) = (V/LY^'^ is the interstring distance for a total length, L, of string in a large 
volume V. We introduced the dimensionless parameters 7,7 and e defined by 7 = l/H^, 
7 = e = 1/H(, where H = R/R is the Hubble parameter; R oc t^^^, with p = 2 in 

the radiation era and p = | in the matter era. With this definition, in the radiation era, 7 
and 7 are identical to the variables used in [Q] and 0, but in the matter era they are half 
as large. In terms of these variables the evolution equations become [jl| 

7 /3 d^i F\ c TGu 

-pt- = -p+ i + -7 + ^e, (3) 

^ 7 V2 2 4/ 2 ' 2 ' ^ ^ 

-ptl = -p+ f «„i + 0-^)-^ + ^^ + TGi2e, (4) 
7 V 2 / wy 2 

-pt- = -p + ftnl H 7^ F 



e 

+A;c7 + TGfiCe. (5) 

The coefficient parameters and functions in @-(^ were fully defined in here we 
recall the definitions briefly. The parameter a = 1 — 2(x^) used in and 0] was written as 
a = ctni+^F, where ani relates to the effect of loop formation and is probably small compared 
to F, which is defined in Eq. (7.17), and is the contribution of stretching. From the 
simulations of Bennett & Bouchet 0], we have, in the radiation- and matter-dominated eras 
respectively, 

ttRD = 0.14 ±0.04 ttMD = 0.26 ± 0.04, (6) 

so F is likely to be almost twice these values. The function G, defined in [|I|, Eq. (7.23), 
also represents a stretching effect; it is larger than F but of the same order, around unity. 
The rate of loop formation is determined by c, which probably lies in the range 0.1-0.5. The 
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effect of gravitational radiation is represented by the parameter FG/x, where G is Newton's 
constant, /i the string tension and F a numerical constant of order 10-100. The probability 
of intercommuting is governed by the constant Xi of order 0.1, and w, defined by the rate at 
which K{1) approaches its asymptotic form for large /, lies between and 1. The function / 
describes the effect of loop formation. It is certainly less than 2c. In it was argued that 
//2c is likely to be only slightly less than unity, but this argument is not correct; in fact 
it could in principle lie anywhere between zero and unity. C is a constant relating to the 
stretching of strings on very small scales, which is probably also small. The two most critical 
parameters are fc, which describes the excess small-scale kinkiness on loops as compared to 
long strings, and C which determines the rate at which gravitational back-reaction smoothes 
the small-scale kinkiness. 

Initially, soon after strings are formed, the gravitational radiation terms in (^-(^ are 
negligible. Then, as we argued in 0], 7 and 7 will reach scaling values, while (if k is not too 
large — see below) e starts to grow. This is the transient scaling regime. The corresponding 
scaling values 7transient and 7transient are found by setting to zero the right hand sides of 
and d). 

This is the only scaling regime accessible to the numerical simulations From 
and (|^) above, assuming that a^i is negligible, we have (the error bars are our estimates): 

C7transient,RD = 1-14 ± 0.04, 

C7transient,MD = 0.26 ± 0.04. (7) 

The second of these is really a purely notional figure, since the transient scaling regime will 
not last into the matter-dominated era. The simulations of |^ also give 

7transient,RD = 7.2 ± 1.4, 

7transient,MD = 2.8 ± 0.7. (8) 

With reasonable estimates for c, this shows that 7transicnt and 7transient are of roughly similar 
magnitude. 

In using Equations (^ and (^, which determine the transient scaling values, we have 
assumed that and the gravitational radiation terms are negligible. There are then four 
remaining parameters: X; ^5 c and J. The latter two are themselves in principle functions of 
7 and 7, given by |ll], Eqs. (6.17) and (6.27), but there are considerable uncertainties in these 
formulae relating to the effects of the small-scale structure, and it seems better therefore at 
this stage to regard them as independent parameters; c at least can be estimated directly 
from the simulations. 

We have explored numerically the solutions of Eqs. (|^) and (|) for the full ranges of 
the four parameters Xj'W^jC and 5=1 — I /2c. Both w and 5 are restricted to lying within 
the range to 1. In general, we find that the value of c is fairly tightly constrained; the 
other parameters less so. In the radiation era, for x = O-l? ranges within which we find 
solutions consistent with Eqs. (|^)-(H) are 

0.10 < c< 0.18, 

0.62 < w< 1, (RD) (9) 
< 5< 0.21. 
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Generally speaking, the smaller values of w or the larger values of c require small values of 
S. In the matter era, conditions are less restrictive; in fact 6 is unconstrained: 



0.08 < c< 0.14, 

0.26 < w< 1, (MD) (10) 
0< 5<1. 

In both cases, lowering the value of x tends to relax the constraints; with x = 0.06, much 
wider ranges of w and 6 are possible — for example, in the radiation era 6 can be as large 
as 0.6. 

The predicted values of 7 are generally towards the lower end of the range allowed by 

(I): 

0.0 <i 7transient, RD < 6.5, 6.5 < 7transient,RD < H, 

2.1 < 7transient,MD < 2.8, 1.7 < 7transient,MD < 3.0. (11) 

It is interesting that 7 tends to be bigger than 7 in the radiation era, but about the same size 
in the matter era. We also note that x ci^nd 7 always appear in (^-(^ in the combination 
X7^; hence lowering x will always increase the scaling value of 7. 

During the transient scaling regime, e grows until the point where the gravitational 
radiation terms become important. The important question is whether then a full scaling 
regime is reached, in which all three of 7, 7, e are constant. If there is such a solution, it is in 
the region of parameter space where e ^ 7fuii ~ 7fuii, which means that the intercommuting 
term in the third equation, xi^ l^i is negligible. If we neglect it, then, assuming we know the 
values of the various constants, it is straightforward to solve for the scaling parameters. As 
we showed in |]l| , the existence of a full scaling solution depends primarily on the magnitudes 
of C and k\ we require C > Ccr > k, where 

p_a,-|(l-4C)F 

and is typically of order one. Unfortunately, we do not know the values of C and k, but we 
shall assume, as seems likely, that these inequalities are satisfied. 

It is important to realize that the simulations performed so far, which neglect gravita- 
tional radiation, may have provided misleading information about the true values of the 
scaling parameters. From the evolution equations, it is clear that when the gravi- 

tational radiation terms become significant, then 7 decreases. Without knowing the values 
of C and k it is not possible to predict the extent of the decrease, but what we can do is to 
examine the relationships between the changing values of 7, 7 and e. 

From Eq. (|^), we might expect that when full scaling is reached the values of 07 and of 
e = FG/ie should be comparable; as e increases, 07 must decrease. The ratio between the 
two is determined by the values of C and fc; by [0, Eq. (9.22), 

e _ TGfie _ Ccr - k 

C7 ~ C7 ~C-Ccr' 
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We have studied the effect of changing e on the values of 7 and 7 predicted by Eqs. (|^) and 
(I). As we noted, as e increases, 07 decreases. It is not so obvious from the equations what 
happens to 7. However, we find that in all cases it too decreases, though by a smaller factor. 

In the radiation era, the dependence of 7 and 7 on e is nearly linear; less so in the matter 
era. For example, if we choose the parameter values x = 0.1, c = 0.1,5 = 0.1, if = 0.8, we 
find the dependence shown in Figure 1. Other parameter values within the allowed range 
give very similar curves. When e = 07, 7 has decreased by nearly a half in the radiation era, 
and in the matter era by about a third. At the same point, 7 has decreased in the radiation 
era by about 25%; in the matter era by 10%. 

From (p!TD, reasonable estimates seem to be 

7full,RD = 3 to 6, C7full,RD = 0.2 to 1, 

7fuii,MD = 1.4 to 2.6, c7Mi,MD = 0.05 to 0.27, (14) 

For observational tests, we need to know how far out we should expect to have to look to 
find a cosmic string. What is most important here is the scaling value of 7 in the recent, 
matter-dominated epoch. With Hq = 2/3t = lOOh km s~^Mpc~^ (0.5 < h < 1), the mean 
interstring distance today would be (7fuii,MD-f^o)~^ = (1200 to 2000)/;.^^ Mpc, with the 
nearest string to us being a distance say half that, i.e. in the range z = 0.2 to 0.5. 

We are also of course interested in the nearest loops, so we turn next to them. The loops 
that are born off the network are usually produced with quite large centre of mass velocities, 
and lose energy both through red-shifting and, at a constant rate, to gravitational radiation, 
Eg^ = — FioopsG/i^, until their rest-mass is reduced to zero. Here Fioops is a constant of order 
50-100. We assume that the fraction of the total energy which will eventually be radiated 
is /, i.e., (1 — /)-Eb is red-shifted away, where i?b is the average energy once it has reached 
a stable non-self-intercommuting form. 

As a first approximation, we imagine that the loop immediately loses the fraction (1—/) of 
its total energy, leaving a static loop of length /i?b//i, which loses its energy by gravitational 
radiation, a fairly good picture as long as the lifetime is long compared to tb- Given this, 
we take the effective length of the loop at birth to be /b = /-f^bZ/i = (k — 1) Fioops Gyutb, 
where we have introduced the dimensionless parameter k. Note that k is identical to the 
1/(3 of Allen and Caldwell and Allen and Shellard At a later time the loop's length 
is / = FioopsG/i(Ktb — t)] it disappears at = Ktb- Numerical simulations indicate k to be 
somewhere between 2 and 10, the latter value certainly satisfying the asumption we have 
made on the role of the initial red shifting. (Even for k ^ 2, the corrections are not large.) 

In a volume V, the rate at which loops are born from long strings is given by. 



uV 

{k — l)FioopsG'/itb 

where 

/C77^ 

^ = 3~ 



(16) 



Taking / to be around 0.7 while 7 and 07 are given by the estimates above, we obtain 

z/RD = 0.2 to 3, i^MD = 0.02 to 0.4. (17) 
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The parameter v is not very well constrained. However, we predict a somewhat lower rate 
of loop formation than for example ref. @]. 

Any possible direct observational test of the cosmic string scenario depends on knowing 
how many long strings or loops we may expect to see within a reasonable distance, and how 
large they would be. We therefore turn to estimates of these quantities. 

The total number density of loops at time t is obtained from A^, allowing for the expansion 
factor: 



where Rif) is the FRW scale factor. Except for the transition region between RD and MD, 
we obtain 



with 



2 K^/^ _ I 

^rd(«:) = 7T ^, ^MD(ft:) = 1. (20) 

This then gives the mean distance between loops, the present distance being 

h-'l-' = ( t = (^^) 90h-' Mpc (21) 



where r2 = Fioops/lOO and fiQ = lO^G/i. Taking umd = 0.1, Fioops = 100, Gfi = 10^^, we 
obtain n^^^^ ~ 200h^^ Mpc. The nearest loop to us is probably at about half this distance. 
If this is the case there are probably a few hundred loops with z < 0.3. 
The total length of all loops in a volume V at time t is given by 



^loops 



^rftb (^:|^^' rioopsG^(«:tb - 1), (22) 



from which we can write the energy density in loops. 



where 



The mean length of a loop is 



Thus at the present time, we have 



Pioops = B{k)^, (23) 



^Rd('^) 3{^3T^ ' 

Bmb{^) =-^ln«:-l. (24) 



1 = ^^rioopsG/it. (25) 
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/ = SMD(/t)r2/i6200/i-^ kpc, (26) 

with a corresponding mass 

Mioop = 1x1 = 5MD(«:)r2/i^4 X 10'^ Mq. (27) 

The median loop size is given by the same expression but with -Bmd replaced by ^^^i- Thus 
we see that the size of a typical loop is larger than a typical galaxy while its mass is that of 
a very large galaxy. In other words, surviving loops are significant objects gravitationally. 

The total gravitational wave energy generated by cosmic strings is the sum of that 
generated by loops and by long strings. In the RD era, we have 

pAt) = f^iAnMi^ + ^isll, t^lr)' (^]\ (28) 

where 



Jtf t \tj \R 



^^is = ^3 • (29) 

Here the lower limit is the time of formation of cosmic strings, or rather the time at 
which they start to move freely and so lose a significant amount of energy to gravitational 
radiation. It is interesting to note that the ratio between the parameters and u which 
govern the contributions of long strings and loops is given by 

^ = (30) 

i.e., essentially the ratio ([T3|). Since 7 decreases as e increases, if is near its upper hmit, 
then u must be relatively small, and vice versa — if more energy is dissipated as gravitational 
radiation less must be available for loop formation. 

The spectrum of gravitational waves emitted by cosmic strings is constrained by the nu- 
cleosynthesis limit on relativistic particle species, as well as the pulsar timing on a stochastic 
gravitational wave background. At the time of nucleosynthesis we obtain 

^gritnuc) = [A^-d{h)v + Z/J lu . (31) 

-J \ 9*,r / if 

Typically, we have T^nc/Ti ~ 10-^^ 9*,nnc/9*,v = 10.75/106.75 ~ 0.1. With these values, 
(pT|) becomes 

f^gr(tnuc) = 1.2 X 10-^6 [^Rd(«;)'^ + ^^Is] ■ (32) 

The nucleosynthesis limit places a limit on the energy density in neutrinos plus any 
relativistic particle species beyond the photons and electrons present at tnuc- It is bounded 
by the equivalent energy density of N^, neutrino species. In terms of the energy density in 
gravitational radiation at nucleosynthesis we obtain 

r^gr < 0.163 X {N^ - 3)aad, (33) 
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where three hght neutrino species have been assumed. According to , the current obser- 



vational hmit on hght (< lOMeV) neutrino species is < 3.1, which yields the bound 

f^gr(tnuc) < 0.02. (34) 

Combining the two expressions we obtain 

^iq[A^d{h)v + z/is] < 15. (35) 

which is not at all restrictive for any values of v and k in the expected range. A more 
stringent limit, Ny < 3.04, has been given recently ||TT|. Using this limit would change the 
bounds in (|^) and (^) to 0.007 and 5, respectively. This bound is still not hard to satisfy, 
but may be somewhat restrictive if either of the terms within the square brackets in ( |53| ) is 
close to its upper limit. 

This though is probably not the tightest gravitational radiation bound. The pulsar 
timing limit requires that the energy density in gravitational radiation, in a logarithmic 
interval at / = (8.2 years)~^ satisfy the inequality ||12 



din/ 

This gives 



< 1.0 X 10"^/i~2. (36) 

-1 



/=(8.2 years)- 



in 



1 /10\ /0.1\3/2 ^1/2 



/Z6<0.2^/— (-)( — ) ^. (37) 



(The details of this derivation will be presented elsewhere.) Note the strong dependence on 
I/, K and h. For h = 0.5 the condition is not particularly restrictive, though it puts bounds 
on u and k, but for = 1.0 it is a very restrictive condition, which could only be satisfied 
by making n and u very small, say k, = 2, u < 0.05. 

The growth of primordial density perturbations generated by strings provides a separate 
constraint on G/i. In the work of Albrecht and Stebbins with strings and hot dark 



matter, the power spectrum of the overdensity fluctuations is approximated by an integral 
in conformal time over the product of the transfer function describing how the initial per- 
turbation spectrum is evolved to the present, and a "form factor", JF, which describes the 
input from the string. It is in this term that the string scaling behaviour is hidden. In 
particular is given in terms of the wavenumber k as 



Briefly, ^ and x correspond roughly to ^ and ^ in our notation, j3 gives the r.m.s. velocity of 
the strings, and E the increase in the surface density of the wakes due to the wiggliness. For 
simplicity, we will assume that the only parameters that are affected are x ^-nd S,- Albrecht 
and Stebbins incorporate the small scale-structure on the network into a renormalized string 
tension /z, and let this value vary depending on the type of string they are attempting to 
mimic. We would like to know in what direction the power spectrum goes if the scaling 
solutions are changed, by increasing x and ^. If xH is doubled (which is roughly the upper 
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limit of what we might expect), then the form factor T is multiphed by 4 and so therefore is 
the power spectrum. Hence in order to fit the large scale structure data, we would naively 
expect that the normalization of [i must go down by a factor of two. For the case where 
they attempt to match the simulations of ref. [Q, they use /J^S ~ 1.2 and (for h = 0.5) 
C/iren = 4.0 X 10~^. Wc would obtain G/iren ~ 2 X 10"^, a value well within current allowed 
limits on Gfi. 

The third observational test which constrains the string parameters arises from aniso- 
tropy observations on the microwave background sky on both large and small angular scales. 
There are few cosmic string calculations of this effect to date, unlike say the situation in 
inflation or global textures. In a recent paper, Perivolaropoulos obtains an analytic 
expression for AT/T in terms of the interstring distance ^. There are many assumptions in 
the calculation, including the use of straight strings and the neglect of loops, which make 
direct comparison somewhat dubious. However, his final result, in our notation, is that 
AT/T oc Gfi'j. Perivolaropoulos quotes a result of Gfi = (2.0 ±0.5) x 10~^. If as we believe, 
7 is somewhat less than the prediction of ref. 0], this should be increased, but only very 
slightly. 

We conclude that the modifications suggested by our analytic study slightly weaken 
the observational constraints on Gfi. For h = 0.5, they are fully compatible with current 
observational limits on this parameter, but for h = 1.0, it is very difficult to accommodate 
the gravitational-wave bound. However, there are still many uncertainties that need to be 
resolved before we can make really definitive predictions. The biggest problem is to find a 
way of estimating the parameters associated with the small-scale structure on the strings, in 
particular C and k which determine the effect on the small scale length of the gravitational 
back-reaction and loop formation, respectively. 
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FIGURES 
FIG. 1. 

Variation of 7 and 7 with e in the full scaling regime (a) in the radiation-dominated and 
(b) in the matter-dominated era. 
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